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ABSTRACT 

The Polytope Escape Problem for continuous linear dynam¬ 
ical systems consists of deciding, given an affine function 
/ : ^ and a convex polytope P C R*^, both with ra¬ 

tional descriptions, whether there exists an initial point xq 
in V such that the trajectory of the unique solution to the 
differential equation 

x{t) = f{x{t)) 
a:(0) = xo 

is entirely contained in V. We show that this problem is re¬ 
ducible in polynomial time to the decision version of linear 
programming with real algebraic coefficients. The latter is a 
special case of the decision problem for the existential the¬ 
ory of real closed fields, which is known to lie between NP 
and PSPACE. Our algorithm makes use of spectral tech¬ 
niques and relies, among others, on tools from Diophantine 
approximation. 

CCS Concepts 

•Theory of computation Timed and hybrid mod¬ 
els; 
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1. INTRODUCTION 

In ambient space R"^, a continuous linear dynamical sys¬ 
tem is a trajectory x{t), where t ranges over the non-negative 
reals, defined by a differential equation x{t) = f{x{t)) in 
which the function / is ajflne or linear. If the initial point 
a;(0) is given, the differential equation uniquely defines the 
entire trajectory. (Linear) dynamical systems have been ex¬ 
tensively studied in Mathematics, Physics, and Engineer¬ 
ing, and more recently have played an increasingly impor¬ 
tant role in Computer Science, notably in the modelling and 
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analysis of cyber-physical systems; a recent and authorita¬ 
tive textbook on the matter is [^. 

In the study of dynamical systems, particularly from the 
perspective of control theory, considerable attention has been 
given to the study of invariant sets, i.e., subsets of R"^ from 
which no trajectory can escape; see, e.g., [inilllElllO]. Our 
focus in the present chapter is on sets with the dual prop¬ 
erty that no trajectory remains trapped. Such sets play a 
key role in analysing liveness properties in cyber-physical 
systems (see, for instance, my- discrete progress is ensured 
by guaranteeing that all trajectories (i.e., from any initial 
starting point) must eventually reach a point at which they 
‘escape’ (temporarily or permanently) the set in question. 

More precisely, given an affine function / : R'^ —^ R"^ and a 
convex polytope P C R"^, both specified using rational coef¬ 
ficients encoded in binary, we consider the Polytope Escape 
Problem which asks whether there is some point xq in V 
for which the corresponding trajectory of the solution to the 
differential equation 

x{t) = f{x{t)) 
a:(0) = Xo 

is entirely contained in P. Our main result is to show that 
this problem is decidable by reducing it in polynomial time 
to the decision version of linear programming with real al¬ 
gebraic coefficients, which itself reduces in polynomial time 
to deciding the truth of a sentence in the first-order theory 
of the reals: a problem whose complexity is known to lie 
between NP and PSPACE [9]. Our algorithm makes use 
of spectral techniques and relies among others on tools from 
Diophantine approximation. 

It is interesting to note that a seemingly closely related 
problem, that of determining whether a given trajectory of 
a linear dynamical system ever hits a given hyperplane (also 
known as the continuous Skolem Problem), is not known to 
be decidable; see, in particular, [aniiiii]. When the target 
is instead taken to be a single point (rather than a hyper¬ 
plane), the corresponding reachability question (known as 
the continuous Orbit Problem) can be decided in polyno¬ 
mial time M- 

2. MATHEMATICAL BACKGROUND 
2.1 Kronecker’s Theorem 

Let T denote the group of complex numbers of modulus 1, 
with multiplication as group operation. Then the function 
0 : R —>■ T given by (l>{x) = exp(27ria:) is a homomorphism 


from the additive group of real numbers to T, with kernel 
the subgroup of integers. 

Recall from the following classical theorem of Kro- 
necker on simultaneous inhomogeneous Diophantine approx¬ 
imation. 

Theorem 1 (Kronecker) . Let 6i,... ,0s be real num¬ 
bers such that the set {6i ,... ,0 s, 1} is linearly independent 
over Q. Then for all tfi,... jtfs £ R and e > 0, there exists 
a positive integer n and integers ni,... ,ns such that 

|n0i — V’l — ^^l| <£,..., |n0s — il)s—ns\ <e . 

We obtain the following simple corollary: 

Corollary 2. Let 9i,...,9s be real numbers such that 
the set {01 ,..., 0s, 1} is linearly independent over Q. Then 

{{4>{n0i ),..., 0(n0s)) : n G N} 

is a dense subset o/T®. 

Proof. Since </> is surjective, an arbitrary element of T® 
can be written in the form ..., cfigfs)) for some real 

numbers tfi,... ,tps. Applying Kronecker’s Theorem, we get 
that for all e > 0, there exists a positive integer n and inte¬ 
gers n\,... ,ns such that 

|n0i — ni\ <£,..., |n0s — ifs — ns\ <e . 

By continuity of (f it follows that ((^(i/ji),. .., (figfs)) is a limit 
point of {{(j>{n9i ),..., (f)(n9s)) ■ n G N}. This establishes the 
result. □ 

2.2 Laurent polynomials 

A multivariate Laurent polynomial is a polynomial in posi¬ 
tive and negative powers of variables zi,... ,Zs with complex 
coefficients. We are interested in Laurent polynomials of the 
special form 

k 

t=i 

where ci,...,cj, G C and ni,i,..., ns,k G Z. We call such g 
self-conjugate Laurent polynomials. Notice that if oi,..., Us G 
T then g{ai ,..., as) is a real number, so we may regard g as 
a function from T® to R. 

Lemma 3. Let g G ..., 2 ^^] be a self-conjugate 

Laurent polynomial that has no constant term. Given real 
numbers 0i,...,0s G R such that 0i,...,0s,l are linearly 
independent over Q, define a function f : R>o —>■ R fey 

fit) = g{(f>{t9i),.. . ,(l){t9s)). 

Then either f is identically zero, or 

liminf f{n) < 0 , 

n—>-oo 

where n ranges over the nonnegative integers. 

Proof. Recall that we may regard g as a function from 
T® to R. Now we consider the function y o (^® : R® —>■ R, 

ixi,...,Xs) gWxi),..., (fixs)). 

We use an averaging argument to establish that either <? o yi® 
is identically zero on R® or there exist x*,... ,Xs G [0,1] such 
that gifjixl),..., 4>ix*s)) < 0. 


Since exp(27rinx)dx = 0 for all non-zero integers n, it 
holds that 

/ ... g{4>ixi),... ,(l){xs))dxi .. .dxs = 0 . 

Jo Jo 

Suppose that g o is not identically zero over R® and 
hence not identically zero over [0,1]®. Then g o (f)‘ cannot 
be nonnegative on [ 0 , 1 ]®, since the integral over a set of 
positive measure of a continuous nonnegative function that 
is not identically zero must be strictly positive. We con¬ 
clude that there must exist (a:i,... , 2 ;*) G [0,1]® such that 
5 f(<^(a;{),..., 0 (x*)) < 0 . 

By assumption, 01 ,... , 0s, 1 are linearly independent over 
Q. By Corollary Q it follows that 

{(0(n0i),..., (finOs)) : n G N} 

is dense in T® and hence has {<j){xi ),..., 4){x1,)) as a limit 
point. Since g o j)’’ is continuous, there are arbitrarily large 
n G N for which 

f(n) = gi(t>in9i),..., (jinOs)) < ..., ^(a:*)) < 0 , 

which proves the result. □ 

Note that this proof could be made constructive by using 
an effective version of Kronecker’s Theorem, as studied in 
[7] and [n], although we do not make use of this fact in the 
present paper. 

We say that a self-conjugate Laurent polynomial g is sim¬ 
ple if it has no constant term and each monomial mentions 
only a single variable. More precisely, g is simple if it can 
be written in the form 

k 

E rij _ —rij 

+ CjZi. , 

where ci,..., Cfe G C, ii,... , 4 G {1,..., s}, and m,... ,nk G 

Z. 

The following consequence of Lemma[3 will be key to prov¬ 
ing decidability of the problem at hand. It is an extension 
of Lemma 4 from [ 6 ]. 

Theorem 4. Let g G ..., 2 ^^] fee a simple self¬ 

conjugate Laurent polynomial and 0i,..., 0s non-zero real 
numbers. Then either 

gi(j>it9i ),..., (j>(t9s)) = 0 for all t gM. 

or 

liminf g{(j){n9 \),..., (j){n9s)) < 0 , 

n—^oo 

where n ranges over the nonnegative integers. 

Proof. Note that if 1,0i,..., 0s are linearly independent 
over Q then the result follows from Lemma Otherwise, 
let { 0 ij,..., 0 ij,} be a maximal subset of { 0 i,..., 0 s} such 
that 1 , 0 ij,..., 0 i^ are linearly independent over Q. 

Then, for some G N and each j, one can write 

N9j = ^m-|-^ni0i,^ , 

where m,ni,... ,nk are integers that depend on j, whilst N 


does not depend on j. It follows that for all j and t G R, 

k 

1=1 

1=1 

In other words, for all j > fe + 1, (f>{N6jt) can be written as 
a product of positive and negative powers of the terms 


One can obtain a closed-form expression for powers of 
block diagonal Jordan matrices, and use this to get a closed- 
form expression for exponential block diagonal Jordan ma¬ 
trices. In fact, if Ji is a fc X fc Jordan block associated with 
some eigenvalue A, then 


/A'* nA"-i (")A"-i 

0 A" nA""^ 


0 0 0 

Vo 0 0 


t n \\n—k+1 
\k-l)^ 

/ n \ \n —fc +2 

U-2;^ 




nA"- 

A" 


1 




It follows that there exists a self-conjugate Laurent polyno¬ 
mial h G C[z^^,... not necessarily simple, but with 

zero constant term, such that for all t G R, 

g{<P{N9it),<t>{Nest)) = h{<t>i9i,t),<p{9i,t)). 

Since l,9i^,... ,9ii, are linearly independent over Q, the re¬ 
sult follows by applying Lemma [3 to /i. □ 


and 


exp(Jit) 


/I 

0 


exp(At) 


t 

1 


0 0 

Vo 0 


Ifc^Tyr 


(fc-2)! 




t 

1 


2.3 Jordan Canonical Forms 

Let A G be a square matrix with rational entries. 

The minimal polynomial of A is the unique monic polyno¬ 
mial m(x) G Q[a:] of least degree such that m(A) = 0. By 
the Cayley-Hamilton Theorem the degree of m is at most the 
dimension d of A. The set o'(A) of eigenvalues is the set of 
roots of m. The index of an eigenvalue A, denoted by v{X), 
is its multiplicity as a root of m. We use v{A) to denote 
max>,go.(A) ^^(A): the maximum index over all eigenvalues of 
A. Given an eigenvalue A G cr(A), we say that u G is a 
generalised eigenvector of A if u G ker(^ — A7)*^, for some 
it G N. 

We denote by Va the subspace of C'* spanned by the set 
of generalised eigenvectors associated with some eigenvalue 
A of A. We denote the subspace of C'* spanned by the set 
of generalised eigenvectors associated with some real eigen¬ 
value by V^. We likewise denote the subspace of C"* spanned 
by the set of generalised eigenvectors associated to eigenval¬ 
ues with non-zero imaginary part by 

It is well known that each vector v £ can be written 
uniquely as u = v\, where ua G Va. It follows that v 

XetT(A) 

can also be uniquely written as u = u’’ J- where G 
and v'^ G VL 

We can write any matrix A G as A — Q~^JQ for 

some invertible matrix Q and block diagonal Jordan matrix 
J = diag( Ji, ... ,Jn), with each block J; having the follow¬ 
ing form: 


In the above, is defined to be 0 when n < j. 

Proposition 5. Let v lie in the generalised eigenspace 
Va for some A G o-{A). Then exp{At)v is a linear com¬ 
bination of terms of the form t" exp(At). 

Proof. Note that, it A = Q~^JQ and J = diag{Ji ,..., Jn) 
is a block diagonal Jordan matrix, then 

exp(j4t) = exp(Jt)Q 

and 

exp(Jt) = diag(exp(Jit),..., exp(Jjvt)). 

The result follows by observing that Qv is zero in every com¬ 
ponent other than those pertaining the block corresponding 
to the eigenspace Va- □ 

In order to compare the asymptotic growth of expressions 
of the form t" exp(At), for A G R and n G No, we define R 
to be the lexicographic order on R x No, that is, 

{g, j) -< (p, m) iff g < p or {g = p and j < m). 

Clearly exp{gt)T = o(exp(pt)t™) as t oo if and only if 

ivJ) -< 

Definition 1. If exp{At)v is not identically zero, the 
maximal (p, m) G R x No (with respect to <) for which there 
is a term t”* exp(At) with 5R(A) = p in the closed-form ex¬ 
pression for exp{At)v is called dominant forb^ exp{At)v. 


/A 1 0 •• • 0\ 

0 A 1 ••• 0 

0 0 0 1 

Vo 0 0 ••• A/ 


Given a rational matrix A, its Jordan Normal Form A = 
Q~^JQ can be computed in polynomial time, as shown in 

0. 

Note that each vector v appearing as a column of the 
matrix Q~^ is a generalised eigenvector of A. We also note 
that the index i^(A) of some eigenvalue A corresponds to the 
dimension of the largest Jordan block associated with it. 


Before we can proceed, we shall need the following auxil¬ 
iary result: 

Proposition 6. Suppose that u G R'* and that 

^ Vx, 

xe<T(Ai) 

where ua G Va. Then vj and v\ are component-wise com¬ 
plex conjugates. 

Proof. We start by observing that 

0 = ^ (uA-uf). (1) 

A6a-(A) 










But if v\ G ker(A —A7)^ then v\ G ker(A —A/)*’, and hence 
vj G Va- Thus each summand v\ — in © lies in Va- 
Since C'* is a direct sum of the generalised eigenspaces of A, 
we must have v\ = vf for all A G o-(A). □ 

We now derive a corollary of Theorem 0. 

Corollary 7. Consider a function of the form h{t) = 
exp(At)u'^, where G V^, with {p,m) G K x No domi¬ 
nant. If h(t) ^ 0, then we have 

—oo < lim inf-— < 0 . 

t-»oo exp(pt)t"* 

Proof. Let 

5R((j(^)) = {rj G 9. : rj i9 G eT{A), for some 0 G R} . 

For each rj G 5R((t(A)) define On = {6 G R>o : rj iO G 
a{A)}. By abuse of notation, we also use On to refer to 
the vector whose coordinates are exactly the members of 
this set, ordered in an increasing way. We note that, due 
to Proposition @ and Proposition [3, the following holds: 

exp{At)v‘' = exp{At) ^ ^ Vn+w + Vn-w 

?7GK((t(A)) 9 € 0 r ] 

= E exp{At)vn+ie 

776SR(cr(A))flee, 

+ exp{At)vn+ie 

= E E C exp(7?t)5r(,,j)(exp(i0^t)) 

77eK(a(A)) j = o 

for some simple self-conjugate Laurent polynomials g(n,j)- 
Note that 

(p,m) = max{(? 7 , j) G R x No : ff(^j)(exp(i6>^t)) ^ 0} . 

The result then follows from Theorem 0 and the fact that 
lim inf--E— = lim inf q(„ ™i(exp(i0„t)). 

□ 

2.4 Computation with Algebraic Numbers 

In this section, we briefly explain how one can represent 
and manipulate algebraic numbers efficiently. 

Any given algebraic number a can be represented as a 
tuple (p,a,£), where p G (!2[a;] is its minimal polynomial, 
a = ai -I- 021 , with 01,02 G Q, is an approximation of a, 
and e G Q is sufficiently small that a is the unique root of 
p within distance e of o. This is referred to as the standard 
or canonical representation of an algebraic number. 

Let / G Z[a;] be a polynomial. The following root-separation 
bound, due to Mignotte [18], can be used to give a value of e 
such that any disk of radius e in the complex plane contains 
at most one root of /. 

Proposition 8. Let f G Z[a:]. If ai and 02 are distinct 
roots of f, then 

I. . I ^ ^ 

|ai 02 ! > ^(d+i)/2jjd-i 

where d and H are respectively the degree and height (max¬ 
imum absolute value of the coefficients) of f. 


It follows that in the canonical representation {p, a, e) of an 
algebraic number a, where p has degree d and height H, we 
may choose ai,a 2 ,e to have bit length polynomial in d and 
log 77. 

Given canonical representations of two algebraic numbers 
a and /3, one can compute canonical representations of a-|-/3, 
ajl, and a/(3, all in polynomial time. More specifically, one 
can: 

• factor an arbitrary polynomial with rational coefficients 
as a product of irreducible polynomials in polynomial 
time using the LLL algorithm, described in HD; 

• compute an approximation of an arbitrary algebraic 
number accurate up to polynomially many bits in poly¬ 
nomial time, due to the work in m-, 

• use the sub-resultant algorithm (see Algorithm 3.3.7 in 
|13p and the two aforementioned procedures to com¬ 
pute canonical representations of sums, differences, mul¬ 
tiplications, and quotient of two canonically represented 
algebraic numbers. 

3. EXISTENTIAL FIRST-ORDER THEORY 
OF THE REALS 

Let X = ( 0 : 1 , ..., Xm) be a list of m real-valued variables, 
and let a{x) be a Boolean combination of atomic predicates 
of the form g{x) ~ 0, where each g{x) is a polynomial with 
integer coefficients in the variables x, and ~ is either > or =. 
Tarski has famously shown that we can decide the truth over 
the field R of sentences of the form if = Qixi ■ ■ ■ QmXmO'(x), 
where Qi is either 3 or V. He did so by showing that this the¬ 
ory admits quantifier elimination (Tarski-Seidenberg Theo¬ 
rem [21] ). The set of all true sentences of such form is called 
the first-order theory of the reals, and the set of all true 
sentences where only existential quantification is allowed is 
called the existential first-order theory of the reals. The 
complexity class 3R is defined as the set of problems hav¬ 
ing a polynomial-time many-one reduction to the existential 
theory of the reals. It was shown in that 3R C PSPACE. 

We also remark that our standard representation of al¬ 
gebraic numbers allows us to write them explicitly in the 
first-order theory of the reals, that is, given a G A, there 
exists a sentence o-{x) such that (j{x) is true if and only if 
X = a. Thus, we allow their use when writing sentences in 
the first-order theory of the reals, for simplicity. 

The decision version of linear programming with canonically- 
defined algebraic coefficients is in 3R, as the emptiness of a 
convex polytope can easily be described by a sentence of the 
form 3xi ■ ■ ■ 3xncr(x). 

Finally, we note that even though the decision version 
of linear programming with rational coefficients is in P, al¬ 
lowing algebraic coefficients makes things more complicated. 
While it has been shown in [I] that this is solvable in time 
polynomial in the size of the problem instance and on the 
degree of the smallest number field containing all algebraic 
numbers in each instance, it turns out that in the problem 
at hand the degree of that extension can be exponential in 
the size of the input. In other words, the splitting field of 
the characteristic polynomial of a matrix can have a degree 
which is exponential in the degree of the characteristic poly¬ 
nomial. 






4. THE POLYTOPE ESCAPE PROBLEM 

The Polytope Escape Problem for continuous linear dy¬ 
namical systems consists of deciding, given an affine func¬ 
tion / : ^ and a convex polytope P C R'*, whether 

there exists an initial point xo £V for which the trajectory 
of the unique solution to the differential equation x{t) = 
f{x{t)),x{0) = xo,t > 0, is entirely contained in V. A 
starting point xo G V is said to be trapped if the trajec¬ 
tory of the corresponding solution is contained in V, and 
eventually trapped if the trajectory of the corresponding so¬ 
lution contains a trapped point. Therefore, the Polytope Es¬ 
cape Problem amounts to deciding whether a trapped point 
exists, which in turn is equivalent to deciding whether an 
eventually trapped point exists. 

The goal of this section is to prove the following result: 

Theorem 9. The Polytope Escape Problem is polynomial¬ 
time reducible to the decision version of linear programming 
with algebraic coefficients. 

A d-dimensional instance of the Polytope Escape Problem 
is a pair {f,V), where / : R*^ —>■ R"^ is an affine function and 
P C R"^ is a convex polytope. In this formulation we assume 
that all numbers involved in the definition of / and V are 
rational^ 

An instance {f,P) of the Polytope Escape Problem is said 
to be homogeneous if / is a linear function and P is a convex 
polytope cone (in particular, cc G P, a > 0 ax G V). 

The restriction of the Polytope Escape Problem to homo¬ 
geneous instances is called the homogeneous Polytope Es¬ 
cape Problem. 


Lemma 10. The Polytope Escape Problem is polynomial¬ 
time reducible to the homogeneous Polytope Escape Problem. 

Proof. Let (/, P) be an instance of the Polytope Escape 
Problem in R'*, and write 

f{x) = Ax -\- a and "P = {a: G R'^ : Bix > bi A B^x > ^ 2 } • 
Now define 

and 



Then {gyV) is a homogeneous instance of the Polytope Es¬ 
cape Problem. 

It is clear that x{t) satisfies the differential equation x{t) = 
f{x{t)) if and only if ( ) satisfies the differential equa- 


tion = g == Q . In general, in any 


^The assumption of rationality is required to justify some 
of our complexity claims (e.g., Jordan Canonical Forms are 
only known to be polynomial-time computable for matri¬ 
ces with rational coordinates). Nevertheless, our procedure 
remains valid in a more general setting, and in fact, the over¬ 
all 3R complexity of our algorithm would not be affected if 
one allowed real algebraic numbers when defining problem 
instances. 


trajectory 


that satisfies this last differential equation, 


the y-component must be constant. 

We claim that (/, P) is a positive instance of the Polytope 
Escape Problem if and only if {g,P') is a positive instance. 
Indeed, if the point xq G R‘* is trapped in (f,P) then the 


point 

Xo 

yo 


is trapped in {g,P'). Conversely, suppose that 


is trapped in {g,P'). Then, since B[ 


Xo 

yo 

-1 


> 0 , we 


must have yo > 0. Scaling, it follows that 

trapped in {g,P'). This implies that yf^xo is trapped in 
(L-P). □ 


We remind the reader that the unique solution of the dif¬ 
ferential equation x(t) = f{x{t)),x{0) = XQ,t > 0, where 
f(x) = Ax, is given by x{t) — exp(At)a:o. In this setting, 
the sets of trapped and eventually trapped points are, re¬ 
spectively: 

T = {xo G R'^ : Vt > 0, exp(At)a:o G P} 

ET = {xo G R”* : > 0, exp(At)a:o G T} 

Note that both T and ET are convex subsets of R'*. 


Lemma 11. The homogeneous Polytope Escape Problem 
is polynomial-time reducible to the decision version of linear 
programming with algebraic coefficients. 

Proof. Let xo = Xq -i- Xo, where Xq € V’" and ajjj G V^. 
We start by showing that if xo lies in the set T of trapped 
points then its component Xo in the real eigenspace lies 
in the set ET of eventually trapped points. Due to the fact 
that the intersection of finitely many convex polytopes is 
still a convex polytope, it suffices to prove this claim for the 
case when P is defined by a single inequality—say P = {x £ 
R"^ : X t> 0}, where > is either > or >. 

We may assume that b^ exp(At)a:Q is not identically zero, 
as in that case 

b^ &cp{At)xo = b^ exp{At)xo 

and our claim holds trivially. Also, if a:o G T, it cannot hold 
that 

b^ exp(At)a:o = 0, 

since b^ exp(At)a;g is negative infinitely often by Corollary^. 

Suppose that xo & T and let {p,m) and {rjyj) be the 
dominant indices for b^ exp{At)xo and b^ exp(At)a:o respec¬ 
tively. Then by Proposition we have 

b^ exp{At)xo = exp(pt)t"*(c 3- o(l)) (2) 

as t —>■ 00 , where c is a non-zero real number. We will show 
that c > 0, from which it follows that Xq G ET. 

It must hold that {ri,j) ^ {p,m). Indeed, if {rj,j) >- 
{p,m), then, as t ^ 00 , 

6" exp(Af)a:o = exp(r,t)f^ , 

a 

but the limit inferior of the termlXlabove is strictly negative 
by Corollary 0, contradicting the fact that xo G T. 





K iVtj) ~ then, as t ^ oo, 

eMAt)xo = eMpt)^" (c + ^ + ®(1)) ’ 

and by invoking Corollary 0 as above, it follows that c > 0. 
Finally, if {r],j) -< (p,m), then, as t —>■ oo. 


Interestingly, the seemingly closely related question of whether 
a given single trajectory of a linear dynamical system re¬ 
mains trapped within a given polytope appears to be con¬ 
siderably more challenging and is not known to be decidable. 

In that instance, it seems that the influence of the complex 
eigenstructure cannot simply be discarded. 


b^ exp(j4t)a;g = exp(pt)t"* • o(l), (3) 

and hence, by m and ®, it follows that 

b^ exp(At)a:o = exp(pt)t"‘ (c -I- o(l)) . 

From the fact that xq £ T and that c 7 ^ 0 we must have 
c > 0 . 

In all cases it holds that c > 0 and hence Xq £ ET. 
Having argued that ET 7^0 iff ET nV^ 7 ^ 0 , we will now 
show that the set ET n is a convex polytope that we can 
efficiently compute. As before, it suffices to prove this claim 
for the case when V = {x £W^ : x [> 0 } (where [> is either 
> or >). 

In what follows, we let [K] denote the set {0, ... ,K — 1}. 
We can write 

b^ exp(At) = exp{rjt)Euf^ j-f , 

(’7,i)6o-(A)x[i/(A)] 

where is a vector of coefficients. 

Note that if cc £ and (p,j) £ {o{A) \ R) x No, then 
j^x = 0 , as j-^x is the coefficient of E exp{rit) in 
b^exp(At)a 3 , and V'" is invariant under exp(At). Moreover, 


Er n V = (H n C) u 


{0} if > is > 
0 if 0 is > 


where 
B = 


Pi {a: £ R'* : uf^j^x — 0} 

( 77 j) 6 (o-(A)\K)x[i/(A)] 


c= U 

(’70')6(o-(A)nK)x [i/(A)] 


{a; £ R"^ : uJ^j^x > 0}n 


P {a: £ R'^ : ufp^^^x = 0} 

(p,m)y(ri,j) 


The set ET n can be seen to be convex from the 
above characterisation. Alternatively, note that ET can be 
shown to be convex from its dehnition and that V'" is con¬ 
vex, therefore so must be their intersection. Thus ET n 
must be a convex polytope whose definition possibly in¬ 
volves canonically-represented real algebraic numbers, and 
the Polytope Escape Problem reduces to testing this poly¬ 
tope for non-emptiness. □ 


5. CONCLUSION 

We have shown that the Polytope Escape Problem for 
continuous-time linear dynamical systems is decidable, and 
in fact, polynomial-time reducible to the decision problem 
for the existential theory of real closed fields. Given an in¬ 
stance of the problem (/, V), with / an affine map, our deci¬ 
sion procedure involves analysing the real eigenstructure of 
the linear operator g(x) := f{x) — /(O). In fact, we showed 
that all complex eigenvalnes could essentially be ignored for 
the pnrposes of deciding this problem. 
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